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y—{ ' Abstract 

o : 

psj ■ The purpose of the paper is to provide a general method for computing hitting dis- 

, tributions of some regular subsets D for Ornstein-Uhlenbeck type operators of the form 

^A + F • V, with F bounded and orthogonal to the boundary of D. As an important 
' application we obtain integral representations of the Poisson kernel for a half-space and 

(f) . balls for hyperbolic Brownian motion and for the classical Ornstein-Uhlenbeck process. 

The method developed in the paper is based on stochastic calculus and on skew product 
representation of multidimensional Brownian motion and yields more complete results as 

■ those based on Feynmann-Kac technique. 

■i— > . 

1 Introduction 

A detailed knowledge of the hitting distribution (equivalently: harmonic measure) of 
a domain for a diffusion with the given generator A is fundamental for solving many 

■ potential-theoretic problems, e.g. Dirichlet boundary problem for a domain or Harnack 
inequality or even boundary Harnack inequality for harmonic functions with respect to 

o : x 

In the paper is we compute hitting distributions of some subsets D for operators of 
the form ^A + F ■ V on subsets of W 1 . It is worth to point out that even in the case 

y—i \ of the classical Ornstein-Uhlenbeck diffusion the explicit formulas for half-spaces or balls 

were obtained only quite recently (see \2\ and |14|). Although the inspiration for our 
work comes from the paper (101 . where the potential theory for bounded sets D and the 

^ , operators 

& '■ 

. 5r ■ 2 

was established, the purpose as well as most of technical tools are here different: instead of 
setting up a general theory, we focus on providing explicit formulas for hitting distributions 
for some important operators of the above type and sets D. The importance of explicit 
formulas is highlited e.g. in the recent papers [7] and [IB], where precise asymptotics for 
Poisson kernel for Bessel diffusions were obtained. Throughout the paper we assume that 
the vector field F in (TjQ) is bounded and orthogonal to the boundary of D. The method 
developed here is based on stochastic calculus and Girsanov's theorem and consists in 



Ia + f-v (1) 
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computing various integral functionals of Brownian motion and representing them in terms 
of special functions. 

We provide a closed formulas for the density function of hitting ditribution, i.e. Poisson 
kernel of a half-space or a ball for the hyperbolic Brownian motion or for the classical 
Ornstein-Uhlenbeck process. The importance of the hyperbolic Brownian motion stems 
from the fact that it is the canonical diffusion on hyperbolic spaces; it has also some 
important applications in the risk theory in financial mathematics, see [11] and |19j . 
Explicit integral representations are crucial in obtaining estimates of Poisson kernel and 
also of the Green function [3], [8]. In these papers, the main tool was the Feynman-Kac 
formula, applied to describe the distribution of a stopped multiplicative functional. The 
present approach, based on methods related to Girsanov's theorem enable us to obtain 
representation formulas for the Poisson kernel, different from those mentioned above. The 
advantage of this approach is seen in Theorem 4, where we obtain the precise asymptotics 
of the Poisson kernel for large values of parameters. Another, worth mentioning result is 
the Theorem 6, where we provide adequate representation of Poisson kernel of a ball. Also 
the formula for the Poisson kernel of a ball for the classical Ornstein-Uhlenbeck diffusion 
is more complete than obtained in |14j (as a series representation only). 

The paper is organized as follows. In Section 2 we provide the general framework for 
the next sections. Throughout the paper we assume that in ([JJ we deal with potential 
vector field F on D, orthogonal to the boundary. Under this assumption, with the aid 
of stochastic calculus and Girsanov's theorem, we establish a general formula for the 
harmonic measure of the set D (Theorem 2). 

In Section 3 we provide a closed formula for the Poisson kernel P(x n ,y) of half-space 
for the hyperbolic Brownian motion on the real hyperbolic space HP (Theorem 3) and 
provide an asymptotic formula for P(x n ,y) (Theorem 4). In Section 4 we provide an 
integral representation of Poisson kernel for centered balls for hyperbolic Brownian motion 
on the ball model D n . We remark here that a similar representation from the paper [5] 
depends on additional conjectures on the zeros of some hyper geometric functions which 
so far remained unsettled. The important tool here, as in the next section, is the skew- 
product representation of the n-dimensional Brownian motion. In Section 5 we provide 
an integral representation for Poisson kernel of a ball for the classical Ornstein-Uhlenbeck 
process (Theorem 7). In Appendix we collect some useful information on Bessel functions, 
hypergeometric and Legendre functions and on skew-product of n-dimensional Brownian 
motion. 

2 Change of measure due to Girsanov theorem 

Notation: For n > 2 we denote by 1™ the n-dimensional Euclidean space, (x, y) denotes 
the standard inner product of x,y G R n and by \x\ we denote the Euclidean length of 
a vector x G BP. The ball with center at zero and the radius r is written as B r = 
{x G IP : \x\ < r}; its boundary, which is the (n — l)-dimensional sphere, is denoted 
by = {x E IBP : \x\ = r} and the spherical measure on is denoted by er™ -1 . 

Furthermore, we write f{x) ~ g(x), when lim x _^b f{x)/g{x) = 1, as x — > b. If for two 
functions / and g there exist constants ci, C2 such that c\ < f(x)/g(x) < C2 for every 
xdDwe will write shortly / ~ g, x E D. 

Throughout the paper D will stand for a domain D C KP with a smooth, connected 
boundary dD and F will be a bounded vector field which is defined on an open Lipschitz 
set U containing D. We assume that F is continuously differentiable up to the boundary 
of U and continuously vanishes on the boundary of U. We further assume that F = W 
on D for a scalar valued function V and call the function V a potential (and F a potential 
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vector field on D). We set F = V = on the complement of U. We say that the vector 
field F is orthogonal to the boundary dD of the set D if for every differentiable curve 
r : [0, 1) dD we have F(T(s)) ■ T'(s) = for every s G [0, 1). 

Lemma 1. Under the above assumptions, if the potential vector field F is orthogonal to 
the boundary of the set D then the potential function V determined, up to a constant, by 
the equation 

W(x) = F{x) 

is constant on the boundary dD. 

Proof. Set xq G dD. The potential V is described by the curve integral 



V(x) = [ F(r)dr + V(x ), 

J<y 



where 7 is an arbitrary continuously differentiable path beginning at xq and ending at 
x. For x G dD we choose 7 to follow the boundary of the set D i.e. 7 : [0, 1] —¥ dD, 
7(0) = xq, 7(1) = x. Then 

V(x) - V(x ) = F F( 7 (s)) • j'(s)ds = . 



JO 

Since dD is connected, we obtain the conclusion. □ 

Throughout the paper we work within the framework of the canonical representation of 
the process, i.e. our basic probability space is the space of all continuous R n -valued func- 
tions defined on [0, 00) with appropriate cr-fields (see [IZ])- The standard n-dimensional 
Brownian motion is denoted by W[t) = (W\{t), . . . , W n (t)). 

Define the process X by the SDE 

dX(t) = dW{t) + F(X(t))dt , (2) 

under the conditions specified above. Then X is a local diffusion on U with the generator 
L = \ A+F(x) • V. Since the field F is bounded, X can be defined as a local semimartingale 
(see, e.g. [El)- Let r be the first exit time of the trajectory from the set D. The harmonic 
measure w x on dD, is defined as the distribution of X{t) under the distribution P x of 
the process X starting at x G D. We define a local martingale M by the formula 

M(t) = f F(W(s)) ■ dW(s). (3) 







Its quadratic variation is then given by the formula 



(M)(t)= \F(W(s))fds. 





We further define the basic object of our study, namely 

N(t) = exp (V(t) - \ (M) (t) \ . (4) 

We now provide the basic formula for the harmonic measure of the process defined by ([2]) 
under some additional conditions. 
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Theorem 2. Under conditions stated above, suppose that X is the process defined by the 
system Assume additionally 

(i) The vector field F is potential and orthogonal to the boundary of the domain D, 

(ii) For every t > 

E x [exp((M) (t A t))] < oo , (5) 

(Hi) 

{N(tAr)}t>o is uniformly integrable. (6) 

Then for x G D the harmonic measure w x has the density function expressed by the 
formula 



w x {dz) = e v(dn)-V(x) E * 



exp (-^ (\VV(W(s))\ 2 + AV(W(s))) ds^j ; W(t) € dz 



where V is the potential function of the field F and V(dD) is its value on the boundary 
3D. 

Proof. According to ([3]) the process 

F{W{t)) ■ dW(t) - ^\F(W(t))\ 2 dt = dM(t) - -d(M) (t) 

is a local semimartingale. Writing, as in @ 

N(t) = exp(M(t) - 1 (M) (i)) 

we obtain, as an application of Ito's formula, that N(t) is a local martingale. If we define 
the measure Q x by 

then, as a consequence of Girsanov's theorem, (W, Q' x ) and (X, P x ) are different descrip- 
tions of the same process, up to time r, see [ID]. Consequently, for a continuous bounded 
function / defined on R n we obtain 

E x f(X(t A r)) = E x (N(t A r); f(W(t A r))) 

Now, the condition ([6]) shows that the expression on the right-hand side converges to 
E x [N(t); f(W(r))], as t — > oo. The left-hand side converges to E x f(X(r)), by the conti- 
nuity of the process X. This indicates that indeed w x has density given by 

w x (dz) = E x [N(t) ; W(t) £ dz] (7) 

We now provide a further description of the function w x . Recall that F is the potential 
of the vector field B. Define 

Z(t) = V(W(t)). 
Applying the ltd formula we see that 

Z(t) - Z(0) = [ VV(W(s)) ■ dW(s) + \ I AV(W(s))ds = 
Jo 2 Jq 

= M(t) - I (M) (t) + ~ (M) (t) + lf AV(W(s))ds. 
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Consequently 



N(t) = exp(M(t) — ~ (M) (t)) 



= exp (z(t) - Z(0) - \f [\VV(W(s))\ 2 + AV(W(s))] ds^j . (8) 

Remark that Y(t) = V(W(r)) = V(dD) since the vector field F is orthogonal to the 
boundary of the set D. Hence, stopping at t A r and taking expectation we get, when 
t — > oo 



exp [|VF(VF(s))| 2 + AV{W{s))) ds^j ; W(t) € dz 



This, together with ([7|), finishes the proof. 

□ 



3 Harmonic measure of hyperbolic horocycle in 

HP 

For every a > we define H a = {x £ W 1 : x n > a}. In this section we consider the 
harmonic measure of the set H a for the operator 

n d 2 d 
&lb = x^y^— -5 - (n-2)x n - — , n ^ 2. 

The motivations for studying this operator comes from hyperbolic geometry. More pre- 
cisely, this operator is the Laplace-Beltrami operator associated with Riemannian metric 
in the half-space model HP of the real n-dimensional hyperbolic space. From geometric 
point of view, the set H a is an interior of the hyperbolic horocycle dH a = {x £ M. n : 
x n = a }- Let (-Bj(i))i=i...n be an n-dimensional Brownian motion on W 1 with the gen- 
erator (and not i.e. the variance E°Bf(t) = 2t. Then the Brownian motion 
on H n , Y = {Yi)i=\,„ni can be described by the following system of stochastic differential 
equations 

f dY x (t) = Y n {t)dB x {t) 

dY 2 (t) = Y n {t)dB 2 {t) , 

{ dY n (t) = Y n {t)dB n (t)-{n-2)Y n (t)dt. 

By the Ito formula one verifies that the generator of the solution of this system is Alb- 
The Laplace-Beltrami operator can be rewritten in the form Alb = 2x n Li, where 

L 1 = ^A + F 1 (x)-V, (10) 
with F\(x) = (0, . . . , 0, (2 — n)/(2x n )). Now, we make change of time. Namely, we write 

A(u)= [ U Y 2 ( S )d S , 
Jo 



and 

If we now write 



a(t) = inf{u > 0; A(u) > t) . 

B k (t)= / Y n (s)dB k (s), k = l,...,n 
Jo 
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then Bk are martingales with mutual variations < B k ,Bi > (t) = 5(k,l)t, k, I = 1, ...n so 
B = (Bfr) is the standard n-dimensional Brownian motion. Substituting 



Y k (t) = Y k (a(t/2)), k = l, 



, n 



we obtain that ([9]) transforms into the following the system of SDE 



dYi(t) = dBi(t) 
dY 2 (t) = dB 2 (t) 

dY n (t) = Z n (t)dB n (t) - (n - 2) 



dt 



2Z n (t) 



Again, by Ito's formula, we see that the operator L\ is the generator of the process 
Y = (Ifc)- Since the operation of change of time does not affect the exit place, the 
harmonic measures of the operators Alb and L\ are the same. 

The potential of the vector field F\ is given by V\{x) = (2 — n) ln(x n )/2. Moreover, it 
is easy to check that the vector field F\ is orthogonal to the boundary of H a . Using ([3]) 
and dH|), we obtain 



M(t) 
N(t) 



n /* dW n (s) 
W n {s) ' 

(n-2)/2 



2 

W n (0) 

W n (t) 



exp 



(M)(t) - 



n 



ds 



ds 



where W(t) = (Wi(t), . . . ,W n (t)) denotes the standard Brownian motion in M. n starting 
from W(0) = x. If we put r a = inf{t > : W(t) g H a } we obtain 



E x exp((M)(t ArJ) 



W n (0) 



E x exp 
E x exp 

(n-2)/2 



(n-2) 2 
4 

(n - 2) 2 



1At„ 



ds 



o 



*Ar a dg 

~n2 



< OO. 



Moreover, we have 

N(t At„) ^ . 

Now, the results of Theorem [2] imply that 

n(n — 2) 



(n-2)/2 



t > 0. 



<{dy) 



n-2 
2 



exp 



ds 



WJsV 



;W(T a )Gdy 



;i2) 



The above-given formula enable us to find the density function P a (x, y), x G -ff a , y G c?il , 
of the measure u x {dy) with respect to the Lebesgue measure on clff a . The scaling prop- 
erty of ?7,-dimensional Brownian motion implies the scaling property for Poisson kernels 
P a (x,y) = a l ~ n Pi(x / a, y/a), x G H a , y G dH a . Moreover, Brownian motion W(t) and 
the set Hi are invariant under translations (x,x n ) — > (x + b, x n + 6), where 6 G M n_1 . 
Consequently, Pi(x,y) = Pi((0,x n ), (y — x, 1)) for every x G Hi and y G dH\. We will 
use these properties in further considerations to simplify the notation. We will denote by 
P(x n ,y) =Pi((0,x n ),(y,l)), where y G M™" 1 and r = n . 

Theorem 3. For every x n > 1 an<i y G M™^ 1 we Ziawe 



P(x n ,y) 



1 



2^-1^+1 



if— l 



t v K v „i(t\y\)dt, (13) 



where v = (n — l)/2. 
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Proof. Observe that the integral appearing in (|12|) as well as the hitting time r depend 
only on the last coordinate of the Brownian motion W(t) = (W(t),W n (t)). Since the 
processes W(t) and W n (t) are independent we obtain 



uf(dy) 



n-2 f 00 



n-2 
2 



where 



^ Xn (ds)=E x 



E°[W(s) G dy]E x " exp 

exp(-|y| 2 /(2a)) 
(2vrs)("- 1 )/2 

n(n — 2) 



n(n — 2) f T ds 



W n {s) 



;t £ ds 



exp 



Vx n {ds) dy, 



ds 



(14) 



WJsV 



yy n 



,T £ ds 



The Laplace transform of fj, x is given by 



E Xn exp 



nm — 2) 



ds 



8 7 W n (s) 2 



w ds 



where q(x) = — n(n — 2)/(8x 2 ) — w. The function ip(x n ) = E Xn e q {r) is a gauge function 
for an appropriate Schrodinger operator based on generator of W n (t). Consequently, (p is 
a bounded solution of the equation 



n(n — 2) 
8x 2 



+ w 



(x) = 0, z > 1 



such that (^(1) = 1. Making substitution y/xip{x^/2w) = <p{x) we reduce the above-given 
equation to 

2x 2 w^"{x\pl~w) + zv^^'^ v 7 ^) - ^ + 2u>£ 2 ^ ij;(xy/2w) = 0, 

which is the modified Bessel equation (|32p with 1/ = . Taking into account the general 
form of solutions of (|32p . the boundary condition and boundedness of (p we arrive at 



K u (x n \/2w) 



w>0. 



(15) 



Since square root can be extended to a holomorphic function on C \ (— 00, 0] and the 
modified Bessel function K v has no zeros in the positive half-plane THew ^ 0, the Laplace 
transform Cfj, Xn (w) can also be extended to analytic function on C \ (— 00, 0]. Moreover, 
using the asymptotic expansion ([33]) we obtain that 



fen Ml < 



-(x n -l)V2^ 1 + E ( x nV^w) 



^2exp (-(x n -l)y/2\ 



w cos ■ 



arg w 



(16) 



for every w G C \ (— 00, 0] such that \w\ is large enough. Note that here argu> G [— tt, tt]. 
In particular, C/j, Xn (w) is bounded for \w\ ^ 1. These properties of C\x Xn and its analytic 
continuation guarantee that we can apply the inverse Laplace transform to (|15p (see |13j . 
Theorem 8.5, p. 267). More precisely, there exist a density function of \i Xn with respect 
to Lebesgue measure on (0, 00) given by the inversion formula 



— : lim 

2iri r-»oo 



l+ir 



C^ Xn (w)e sw dw. 



l—ir 
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To compute the limit we integrate the function f s (w) = £fi Xn (w)e sw over rectangular 
contour surrounding the branch-cut of f s which is negative real axis. Let T be the 
positively oriented contour consisting of four horizontal segments 71 = [— r + i/r,i/r], 
72 = [— r — i/r, —i/r], 73 = [— r + ir, 1 + ir], 74 = [— r — ir, 1 — ir], three vertical segments 
75 = [— r + i/r, — r + ir], 76 = [— r — i/r, —r — ir], 77 = [1 — ir, 1 + ir] and a semi-circle 
78 = {\w\ = l/r,Rew > 0}. The formula (f3l|) implies that f s is bounded for small w 
such that Rew > 0. Consequently, the integral over 73 tends to zero when r — > 00. The 
boundedness of £[i Xn (w) for large w implies that for r ^ 1 and every s > we have 



+ / ) f s (w)dw 

75 ^76 



< 2 sup Cfi Xn (w)re rs -)• 0, 
IH>1 



as r — )• 00. Finally, using (fT6|) . we obtain 



( f + [ ) < 4exp(-(x n - 1) cos(37r/8)V^) / 

VJ73 ^74/ ■/— 1 



e _su du -> 



+ / ) / fl Hdu; 

'73 ^74 

as r — )• 00. The Cauchy Theorem together with previous considerations and (|35p give 

Hx n (s) = ^— lim [ Cfi Xn {w)e sw dw = ^— lim ( [ + [ \ f s {w)dw 

K u (-iV2tx n ) K v {i^/2~tx n ) 



2vri 







K v {-iy/2t) K v (i\f2i) 



e~ st dt 



K„(iV2t Zn ) 



K u {iV2t) 



Im " v ) te~ st ' z dt 

7T ./ V Ai/^t) 



From HSU) and (J3TJ) is easy to see that J^(t) + Y?(t) ~ t" 1 V t~ 2l/ and 



<C(l + i), t>0 



for some constant C = C(x n ) > 0. With the use of (|36|) we verify 

\J v (t)Y u (tx n ) - J v (tx n )Y v (t)\, ( f°° e-\y\ 2 l^e- st2 ' 2 ds 



dt 



roc 

(t + l)t v K v „i(t\y\)dt. 

Jo 



The last integral is finite by (|33j) and (|M|) . Consequently, by the Fubini's theorem we 
obtain 

p , v < Z 00 J,(i)y,(to n )-J,(to n )y y (t) f / z" 00 e-\y\ 2 i^e- st2 ' 2 d s \ 



1 



This ends the proof. 



J v {t)Y v {tx n ) - J„(tx n )Y u (t) 
J 2 »(t)+Y?(t) 



t u K u ^(t\y\)dt. 



□ 



S 



The integral formula presented in Theorem can be used to obtain the asymptotics 
of the Poisson kernel P(x n , y) as well as its sharp bounds for small x n and large \y\. Note 
that results given in the next theorem cover those obtained in [5] (see Theorem 4.9 and 
Theorem 4.10, compare also with Theorem 5.3 in 0). Moreover, the formula (|13p allows 
to omit very laborious and sophisticated computations used in [4] to examine the behavior 
of P(x n ,y) when \y\ tends to infinity. Our approach is simpler and gives more general 
results. 

Theorem 4. For every xq ^ 1 we have 

. r(n/2) ^— 1 k x n — 1 , , 

P ( X n,V) ~ 2 n- 27T n/2 Z^ X 0- y2^=2' X n ^ X , M ^ OC. (17) 

fc=0 |y| 

Moreover, for every yo > we have 

P ( X ™'^ ~ 7r (n + 3V2 / 2 (n-5)/2 ' ~ ^ *n 1, Ifo, (18) 

where 

s u K„_ 1 (sy )ds 



c(yo) = bo 



l-n 
2 



o J2 (s) + y i 2 (s) - 



Proof. Making a substitution t\y\ = s in (|13p we can rewrite the Poisson kernel in the 
following way 



P(x n ,y)\y\ 4l/ = 

X n ~ 1 ~ 



where 



a(xt) - 1 JuWvW - Mtx)Y v (t) x>lt>0 

Since s 3u K u _i(s) is integrable on [0, oo) and ([40]) gives boundedness of 
x n < R, we can apply Lebesgue dominated convergence theorem to get 



9v [X n , |?,| 



for 



r P{x n ,y)\y\ Av 1 r aw f ^ 2r(2i/)r(i/ + i) 



and 

P(x w ,2/)M^ x r°° n{xf - l)s^ K v ^{s)ds 

(x,| w |)™* 0l oo) x n -l 7r«^l2"-iy o 2 2 "(x -l)r(i/)r(i/ + l)sg 

7T-(X " l)r(z.)x^" 1 TT"^-^^) 

whenever xo > 1. Here we used formula (|38p from Lemma [8] and relation (|37|). The 
duplication formula for gamma function gives (|17j) . In the same way, using (|38p . we get 

lim P(xn,y)\y\ Av _ 1 f x 2s v \y \ 2v K v ^{s)ds 



1 ^ v lyol; ^ ^ V lift) I 



which, by substitution s = tyo, proves JT8J). □ 
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As a consequence of Theorem [J] we obtain the following sharp bounds for the Poisson 
kernel for small x n and large |y|. Similar results have been obtained recently (see Theorem 
11 in [7]). The results are more general (there is no restriction for x n and |y|) however 
the methods of proof are much more complicated. 

Corollary 5. We have 

X — 1 

P(x n ,y) « |J| 2n -2 ' Kx n ^2,\y\^l. 

Proof. Existence and positivity of the limits proved in Theorem U] imply that for every 
x £ [1,2] there exist e x > and Y x > and strictly positive constants c\{x),C2{x) such 
that 

x' — 1 x' — 1 

for every x' 6 [1,2] and y £ R n_1 satisfying \x — x'\ < e x and \y\ > Y x . Since the 
family of intervals {x — e x ,x + e x ) is an open cover of [1,2], we can choose the finite 
subcover {{xi - e Xi ,Xi + e x .),i = l,...m}. Putting Y = max{Y Xi : i = l,...,m}, 
ci = min{ci(xj) : i = 1, . . . , m} and c 2 = max{c2(xj) : i = 1, . . . , m} we get 

Cl ^|2^2 <^n,y) ^ C 2^2^2 ( 19 ) 

for every x n G [1, 2] and \y\ > Y. Observe that formula (JT3J) implies positivity of P(x n , y). 
Moreover, by (fT8|) . we get that the function P(x n ,y)\y\ 2n ~ 2 /(x n — 1) can be continuously 
extended to the strictly positive function on the compact set [1, 2] x [1,Y] and consequently 
(|19|) is true also for |y| G [1,^] (with possibly different constants c\ and a?). This ends 
the proof. □ 



4 Harmonic measure of hyperbolic balls in W 1 

In this section we consider the harmonic measure of the balls associated with the operator 

i=l ' i=l 

This operator appears naturally as the Laplace-Beltrami operator on the ball model of 
the real hyperbolic space B" (see [5] for more details). 

In particular, if B = (Bk) is the standard n-dimensional Brownian motion then the 
following system of SDE: 

dY k (t) _ dBk ^ + 2 ( n _ 2 )Y k (t)dt, k = 1, . . . , n. 



l-|F(t)| 2 



describes a diffusion with values in the real hyperbolic space B n , with the generator 
2 Alb- As in the case of half-space model, we perform a change of time defined by: 

pu 

A(u)= (1 - \X(s)\ 2 )ds, and a(t) = inf {u; A{u) > t} . 
Jo 

Then the process defined by 

Y k (t) = Y k (a(t)) 

is the diffusion with values in B n with the generator 
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As in the case of HP in Section 4, the harmonic measures of the operators Alb and 
L 2 are the same. 

We consider now the harmonic measure of a ball B r = {x E R n : \x\ < r}, r < 1, 
supported on the boundary of B r which is the sphere S"™" 1 of radius r. We denote by 
P r (x,y) the Poisson kernel of B r , i.e. the density of the measure u)® with respect to the 
(n — l)-dimensional spherical measure cr™ -1 . As in the previous section we can write the 
Laplace-Beltrami operator in the form i — U ; L 2 , where 



Lo 



1 



A + F 2 (x) ■ V, 



|2\2 



/2 implies that the harmonic functions on the 



with F 2 (x) = yq^rfXi, . . .,x n ). 

The positivity of the factor (1 - 
ball B r for the operators (|20p and L 2 are exactly the same and consequently the harmonic 
measures coincide. 

Moreover, the vector field F 2 is orthogonal to the sphere 5 1 ™ -1 and its potential function 
is V 2 {x) = ln(l - \x\ 2 ). We denote by r r = inf{t > 0; W(t) B r } the first exit time 
of Brownian motion W(t) from a ball B r . Note that r r depends only on the Euclidean 
norm of W. 

The martingale M, related to the vector field F 2 , and its quadratic variation are 



M(t) 



{n -2) 



'^'W (M)(t) = (n-2f 



l-\W(s)\* 

Observe that the condition ([5]) is fulfilled in this case, since 



\W(s)\ 2 ds 

(f^WW?' 



E^exp 



(n 



tATr \W(s)\ 2 ds 
(l-\W(s)\ 2 ) 2 

By (|8j), the kernel iV(i) is of the form 
'1-|^(0)| 2 



N(t) 



l-\W(t)\< 



n — 2 
2 



exp 



< E x exp 



n(n — 2) 



(n-2) : 



,r 2 (tAT r ) 
(1-r 2 ) 2 



< 00 . 



ds 



(1-|^(,)| 2 ) 2 



and it is now evident that {N(t A r r )}j>o is uniformly bounded in t so the condition 
holds. Applying Theorem [2] we obtain 



Wr(dy) 



n-2 

2 



exp 



n in 



ds 



[l-\W(s)\ 



2\2 



\W{r r ) E dy 



(21) 



From now on we assume that i/0. For x = 0, from the rotational invariance of the 
Laplace-Beltrami operator, we easily obtain that uj x is just ct" _1 /(T™ _1 (S'™ _1 ). Recall the 
skew-product representation of the Brownian motion 

W(t) = flM(t)9(AM(t)), 

where BS- V > is the Bessel process with index v = ¥■ — 1 staring from \x\ and is spherical 
Brownian motion on S"" -1 independent from bS u ' (see Appendix). Using the fact that r r 
depends only on BS V / we get that W(r r ) = R^ u \T r )@(A^ u ' (r r )), where 6 is independent 
from R( u \r r ) and (r r ). Applying this decomposition to formula (|2ip we get 



n-2 
2 



PFI(G t G dy)! 1 



exp 



q(RM)ds) ]AM(r r ) G dt 
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where q(y) = — 2 "i_y^ p • Rotational invariance of spherical Brownian motion implies that 
the harmonic measure is axially symmetric with axis x. As a consequence we get that 
its density P r {x,y) depends only on the cosine of angle between starting point x and the 
point y, i.e. 

P r (x,y)=P r (x, {X ' V) 



\ x \\y\ 

If we consider the sets of the form A = {ij G S™^ 1 : G (a, 6)}, where —1 < a < b < 1 
we get (for a definition of the process S see Appendix) 

P^(e t eA) = P 1 (S t G (a, 6)) = f b p?(l,z)m(dz) =2 I pf (1, z)(l - z^'^dz, 

J a J a 

where is defined in (|50p . From the other side, using the spherical coordinates we obtain 
<4(rA) = ! P r {x,y)da 1 }~ 1 {y) = r n - 1 (j 1 }- 2 {S™- 2 ) I P r (x, cos <f>) sin n " 2 d<j) 

J rA J cos (j>£(a,b) 

P r {x,z)(l _ z 2 )(™-3)/2 ^ 



mr™* r n 1 '''' 



r(B±i) 

Comparing both sides we get the following formula for the Poisson kernel P r (x,y) 



vr—nr™- 1 V 1 - r / Jo \ \ x \\y 

where 



/x y (dt) = E y 



exp 



(i? M (s))ds J ]A {u \t t ) G 



, y€(0,r] 



The formula for pf can be computed from the appropriate formula for the transition 
density function for 0, which is given in terms of spherical harmonics and that approach 
leads to the series representation for P r (x,y) presented in [5]. However, we want to 
compute the Laplace transform of pf which is so called A-Green function of the process S 



G x (x, 1) = / e- xt pf (1, x)dt, x G (-1, 1), 
J o 



and we do it directly. From the general theory (see for example |3] Chapter II for short 
resume) a function G\ is described by solutions of the second-order differential equation 



1 — x 2 it/ \ n — l 



u"{x) xu'(x) = Xu(x), a; €(-1,1). (23) 

z 2 

Note that the expression on the left-hand side is just Qu(x), where Q is just the generator 
of S described in (|49|) . More precisely, we have 

y A (l)V> A (x) 

G\{x,l) = , a; €(-1,1), 

w x 

where ip\ is a decreasing and tp\ is an increasing solution of (f23j) such that </?^~(l) = 
and tp^(-l) = 0. The boundary conditions for the derivatives follows from the fact that 
non-singular points —1 and 1 are reflecting. Here f + and /~ denote the right and left 
derivative with respect to the speed function s(x). The Wronskian w\ is given by 

w \ = i>l{x)tp\{x) - ?jj\(x)ip^ (x) 
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and it does not depend on x. Putting x = 1 in the above-given formula and using the 
boundary conditions we obtain that 

A( ' j Va(i)va(i)- ^(i)^(i) KO)' ( j ' 



This implies that G\(x, 1) is uniquely described as a solution of ([23]) such that u + (— 1) = 
and u~(l) = 1. Making a substitution u(x) = f (z) with z = in the equation (j23|) we 
reduce it to the following hypergeometric equation 

z(l - z)f\z) + - (n - /'(*) - 2A/(z) = 0. 



with a = ^-^(A), /3 = 2f2+4(A), 7 = ^ , where A(A) = \ ^/(n - 2) 2 - 8A. Consider 
the function 

The above-given computation implies that the function h\ is a solution of (|23p . Using 
(|43|) and (|44|) we compute the derivative of this function with respect to the scale function 
s(x) in the following way 



n-l 



d , , . . 9, n-i 2A _ / n ,,,x n .... n+1 l + i 



l-x 2 ) — —h x (x) = (l- x ^ — —- 2 F 1 (--A(X),- + A(\y, 



dx K ' v ' n-l \2 v " 2 v " 2 ' 2 

n+1 . 

2 2 A., .n-i ,_, /l .... 1 .... n + 1 1 + x 
= —( l + X ^^{2 +A ^2- A ^—- — 

The first equality and the fact that the hypergeometric function 2 F\ is equal to 1 at zero 
implies h~t(—l) = 0. Using the second equality and (|42p we obtain 



r(a=2_A(A))r(2f? + A(A))' 



Moreover, using once again (j44J) and the definition (|46[) we can express function h\ in 
terms of the Legendre function of the first kind 

3 — n 

, , \ fl — x\ 2 ^ (\ 1 n-l 1 + x 

M*) = (— J 2F 1 f < --A(X),- + A(Xy,— ; — 

(l-x 2 )^ /n-l\ 2^ 
Finally we have just obtained that 

3—n 3—n 

G A 0z,l) = B„(A)(l-x 2 ) — P A( 2 A) _i(-x), 

where 

„ m r (ag - .4(A)) r (ag + .4(A)) 
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The second part of the formula ([22]) relates to the measure \i y . Observe that /j, y 
depends only on the Bessel process which is one-dimensional diffusion. For every 

w ^ the Laplace transform Cfi y (w) is 



E^exp 



n(n 



ds 



(l-(/2M( S ))2)5 



w 



ds 



o (i?M(s)) 5 



We g (T r ), 



~~ ^he f unc ti° n ^(y) = E y e s (T r ) is by definition the gauge 



where g(y) = 

function for the Schrodinger operator based on the generator of the process flM and the 
non-positive potential g. From the Feynman-Kac formula ip is a solution of the Schrodinger 
equation. Using (|48[) we obtain that <p is a bounded solution to the following second-order 
differential equation 



2^ + ~Ty 



TliTl — 2 J TV 

2(1-^)2 +^)^//) -^ //£[()./■), (2.1) 



with the boundary condition f(r) = 1. Substituting y(y) = y n ' ip yjZy^j we obtain 

1 + 2/ 2 ' 



-n/2 



■0 



+ 



4y 2-n/2 /l + y 2 



(1-2/ 



2\2 



n /n 



fey 3 -"/ 2 .„ n + y' 

+ ~ —IP 



(1-2/ 



1 - y 2 ^ 

2 



i-r 



+ 



4 l-n/2 / 1 , 

+ 7T3(3-«+(Hl)|/V / ' 



i - r 



Putting the above given formulas to the differential equation 
by the factor - T^z^p give 



and dividing both sides 







(1-2/ 2 ) 2 



2(i-y 4 )y A + 2/ 2 



1-1/ 



2^2 



+ 



(i-y 

4y2 



Moreover, putting z = jzJjp and using the equality 1 
following differential equation for i/j 



2\2 



n 



+ 2w 



n(n-2)\ , /l + y 2 



l-y2 

= / 1 ~ 4^ a^a lead to the 
(l-irr 



(1 - z 2 W(z) - 2z^\z) +(v(u+ 1) 



A(-wf 



if)(z) = 0, z > 1. 



with i/ = § — 1 and A(— w) = \ y/{n — 2) 2 + 8w. This is the Legendre's differential 
equation (j4"5j) . Thus, the general solution of (|24"|) is given by 



V 1 - 2/ 



1 + 3/' 



2/ G [0, r] 



where c\ and C2 are absolute constants. Using (|4T)j) and (|47p one can easily check that 
the function y~ u P u A ^ w ^ ( \ +v A ] is bounded on the interval 



1, jz^z ) in contrast to the 



function y U Q U A ^ w ^ ( jzjfp j, which is unbounded in the neighborhood of 1. Thus C2 = 
and the boundary condition ip(r) = 1 gives 



\ n/2-1 P- A ^ f 1±|*£ 



Ixl < r,w ^ 0. 



(25) 
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Now observe that for every complex number w such that Re(w) > — \ 



\L^ A iw)\ < EN 



/ n(n-2) f Tr ds , . f Tr ds 

exp / 7-7— - Re(™) 



2 /o (1 - (i?^) 2 ) 2 (i?^) 2 , 
(n-2) 2 p ds \ M (v 2 r* ds \ fry 



< ENexp — 7-— = E |x| exp -7- 



/o (i?^) 2 ; v 2 ^° {r^yj Kx 

The last equality follows from (see [3] 2.20.4 p. 407) 

In particular £//| x .|(— v 2 /2) is finite. This implies that the formula 



oc 



C^w) = / e- wt fi lxl (dt) 
Jo 

defines a holomorphic function in the complex half-plane Re(iu) > —v 2 /2. Moreover, for 
\z\ < 1 the function 2-^1(0, /?; 7; z)/T(a) as a function of a is analytic function on C. 
Using this fact and the representation of Py in terms of hypergeometric function 2-P1 we 
deduce that the function on the right-hand side of (|25p is a meromorphic function in the 
half-plane Re(io) > —u 2 /2. In fact, the equality (|25[) implies that the ratio of Legendre 
functions is analytic for Re(io) > —v 2 j2 and consequently the function in the denominator 
has no zeros in this region. Compare this result with Conjecture 5.2 in [5]. Moreover, we 
have just proved that (|25p holds whenever Re(tu) > —v 2 /2. 
Now let c = — 16 . We have 

roo _ / /_ „,\\ r oo _ / /_ „,\\ / -1 /-c+ioo 



dtdz. 

Taking into account the previously found formulas for the Laplace transform Cfi\ x \{z) and 
the Green function G z (l,x) we finally obtain 

Theorem 6. For every x G B r , x/O and y G S*™ -1 i/ie Poisson kernel P r (x,y) is given 
by the following formula 



2T(*P) fl-\x\ 2 r Vsin 2 ^ „ ,.,p- 



1 1 'iiii t x B n (z)P .J . , (— cos o?)dz, 

w/iere A(z) = ±^(n-2) 2 -8*, B n (z) = r (^~ffi)) r (g +A < 2 )) ? c = _^2/ 4 fflnd ^ zs fln 
angle between x and y, 



5 Harmonic measure of Ornstein-Uhlenbeck pro- 
cess 

As in the previous chapter for fixed r > we denote B r = {x G K" : \x\ < r}. Consider 
a vector field Fz(x) = Xx, where A > 0. F3 is a potential vector field with the potential 
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function Vs(x) = X\x\ 2 /2 and as in the previous cases it is orthogonal to the boundary of 
B r . The corresponding martingale M and its quadratic variation are 



M W = A E / Wi(s)dWi(s), (M) (t) = / X 2 \W( 
i=1 Jo Jo 



s)\ 2 ds. 



The validity of ([5]) in this case follows from 

rlAT r 



E x [exp{M) (t/\r r 



E 2 



exp 



X 2 \W(s)\ 2 ds 



o 



where r r = inf{t > : W(t) G' B r }. Since 



^ E x [exp[X 2 r 2 (t AtJ < oo, 



(X 2 \W(s)\ 2 + 2nX)ds 



all the assumptions of Theorem [2] are satisfied and consequently we obtain that the har- 
monic measure w x .(dy) of B r for the operator 



L 3 = ^A + Xx-V 



is given by 
w x r {dy) = exp A ( r2 ~N 2 V 



exp 



(A 2 |W(s)| 2 + 2nA)(is ; W(T r ) € dy 



Computations in this case mimic those introduced in the previous section so we omit 
some details and present only a sketch of the argumentation. For i/O the skew-product 
representation of the Brownian motion allows us to write 



(dy) = et( r2 -N 2 ) / PT=T(e t €dy)ll*l 



exp 



G (it 



where g(y) 



A 2 UI2 



x| — nA and consequently the Poisson kernel P r (x,y) is given by 



where 



P r (x,y) 



Vy{dt) = E y 



2r(2±i) 



(n-l)/2 nr n-l 



e |(r- 2 -k| 2 ) 



7T 



exp 



o (R^) 2 



G 



, y€(0,r]. 



As previously, the Laplace transform C/j, y (w) given by 

/ A 2 



E y exp 



-(i?^) 2 -nA- 
2 (R { ? ] ) 2 , 



ds = E%(r, 



g\'rj, 



(26) 



where = —^-y 2 — nX—-^, can be identified (by applying the Feynman-Kac formula) 
as a bounded solution of the following Schrodinger equation 



y 2 + nX + ^ ) v :(//) - (). ye (()./■). > U. 
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with the boundary condition <p(r) = 1. Setting <p(y) = y n ^ 2 f(Xy 2 ) we reduce this 
equation to the Whittaker equation 



/"(*) + /(*) 



1 n f n(n — 4) ui\ 1 



4 2x V 16 2 / x 2 



0, 



with parameters k = — ^ and ^, = V(" 2 ) Consequently 

¥>(y) = y~ f [ciM(fc, /z, Ay 2 ) + c 2 W(/c, ^, Ay 2 )] , 

where M and W are Whittaker functions (see pQ 13.1.32, 13.1.33 p. 505). The boundedness 
of (/? implies that c 2 = and the boundary condition <p(r) = 1 gives 

— n 

M^ArT (27) 

If we look at ([26]) . the probabilistic definition of C/j, y (w), by the same argument as pre- 
viously gives that Cfi y (w) can be extended to holomorphic function in the complex half- 
plane Re('w) > Since the Whittaker functions are well defined in this region the 
formula (J27J) is also satisfied in this region. As before we use the Laplace inverse formula 
and for c = — ^ n 7? obtain 

16 



poo 


Pt ( 1, 


Jo 




1 


p c+ioo 


2vH 


J c—ioo 


1 


/•c+ioo 


2tH 


J c—ioo 



- r +lCC 'cfi ]x] (z)G t dtdz - 
KiJc-ioo V \ x \\y\J 



Theorem 7. For every x 6 -B r , a; 7^ and y 6 S 1 " 1 t/ie Poisson kernel P r (x,y) is given 
by the following formula 



n-1 

7r 2 nr 



where (p is an angle between x and y, A(z) = \y/{n — 2) 2 — 8z and 

_ r(^-A(z))r(^ + A(z)) 

On\ z ) — n+ i n 

2— ^r(^fi) 

6 Appendix 

For convenience of the Reader we collect here basic information about Bessel functions, 
hypergeometric functions and other special functions appearing throughout the paper. 
Mainly we follow the exposition given in [I] and [12], where we refer the Reader for more 
details (see also [TS] and [T5j). 

6.1 Bessel functions 

The Bessel functions J v {z) and Y v (z) are independent solutions of the Bessel equation 

z 2 y"(z) + zy'(z) + (z 2 - v 2 )y(z) =0, v € R. 
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The Wronskian of the pair (J u (z),Y u (z)) is equal to 2/(irz) (see [18J, p. 113). The deriva- 
tives of Bessel functions can be expressed by itself in the following way 

J' u [x) = J v - X {x)--J v (x), x>0, (28) 
x 

Yl(x) = n_i(x) - -Y v (x), x>0. (29) 

x 

For every i/>0we have (see [18] 5.16 p. 134-135). 

r v 2 v T(v\ 1 

J„(z) —, r , YJx) — — , X^0+. (30) 

v 1 2T(i/ + l) v ; vr s"' v ; 

/ 2 w i / 2 i 

M x )~\ — cos(x- — --) , y Q (x ~ \ — sini- — - -), x -)-oo.(31) 

V X7T 2 4 V X7T 2 4 

The modified Bessel functions are independent solutions to the modified Bessel equation 

z 2 y"(z) + zy'(z)-(v 2 + z 2 )y(z)=0. (32) 
The following asymptotic expansion for the function K v holds (see [15] 8.451 (6)) 



KJz) = J—e- z (l + E(z)), \E(z)\ = Oflzr 1 ) as \z\ -> oo. 
V 2z 



(33) 



whenever | argz| 37r/2 . The behavior of K v near zero is described by (see [I], 9.6.9) 

Ku{z)~ ^- L , Rez>0. (34) 

The connection between modified Bessel function of purely imaginary argument and Bessel 
functions is given by 

K„(ix) = -ye-'"* /2 (Mx) - iY u (x)), x > 0. (35) 
Finally, we recall the integral representation of K v (see [15], 8.432 (7)) 

Kt(z) = Z — J™ exp (~ t+ f^ t-^dt, z>0,$£R, (36) 
as well as the formula ([IS] 6.561 formula 16, p. 676) 

j~ x»K v {x)dx = 2^r r {^f^) , M + 1 > - > 0. (37) 

For every ^ > 1/2 we introduce the following function of two variables 

a( x t)- 1 Mt)Y v (tx) - Mtx)Y v {t) x>lt>Q 

It is obvious that g v is a continuous function on (l,oo) x R + . However, the most crucial 
for the considerations given in Chapter [3] are the following asymptotic properties of g v . 

Lemma 8. Set x$ 1, to > 0. Then 

lim gJx,t) = = J"Vw 0X ° \ z , (38) 

2 
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Proof. From the Lagrange theorem, there exist 9i,9 2 £ (^, x ), depending on x and such 
that 

J v (t)Y v (tx) - J v (tx)Y v (t) _ j M Y v (tx) - Y v {t) J v (tx) - J„(t) 

tx — t tx — t tx — t 

= J v {t)Yl(t6i) - J' v {te 2 )Y u {t). 

Obviously, when x — > 1 then 6\, 9 2 also tend to 1. Furthermore, since the Wronskian of 
(J u (z),Y u (z)) is 2/nz we get that 

r Mt)Yu(tx) - J u (tx)Y v (t) , , 2 

(x,t)-*(l t t ) IX — t TTlQ 

which proves (|39p. If we use the recurrent formulas for the Bessel function derivatives 
and ([291) we obtain that J u (t)Y/,(t9i) - J' u {t9 2 )Y u {t) is equal to 

Multiplying the last expression by t, taking (x,t) — > (1,0) and using ([30]) it is easy to see 
that first two summands tends to zero and the last one tends to 2 /it. Since, by (j30|) . we 
have lim^ + t 2 "(J 2 (t) + Y 2 (i)) = 2 2v Y 2 (u)/t: 2 and consequently 

which is (|55|) for xq = 1. For xq > 1 relation (|55|) follows directly from (|5U|) 

lim (j t) = nx ° ~ X ° ) = n ^k=o x o 

(*,*) =Kx .o) ' ' v2*»T*(v) 2 2 "T(v)T(v + l)xg ' 

□ 

Note that Lemma [8] together with (|30[) and (|3ip imply that the function g v can be 
extended to the continuous function on [l,oo) x [0,oo) which is bounded whenever x is 
bounded, i.e. for every R > 1 there exists C(R) > such that 

\g v (x, t)\ < C(R), (x, t) G [1, R] x [0, oo). (40) 

6.2 Hypergeometric Function and Legendre Functions 

For 7 ^ -1,-2,... the hypergeometric function is defined by 



2 F l(a ,ft 7 ;,)=X;^§4*, N <1 



fc=0 

Here (a)& = T(a + k)/T(a). Function 2-2*1 is a solution of the hypergeometric equation 

z(l - z)u"{z) + [7 - (a + + 1)*K(*) - a/3u(z) = (41) 
regular at z = 0. Whenever Re (7 — a — /3) > we have (see [12] vol.1 p. 104 2.8(46)) 

„, fl u r(7)r(7 - a - (3) 

2 Fi(a,P;j 1) = — — — . 42 

r 7-q r 7-^ 
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The derivative of 2 F X is given by (see [T2] vol.1 p. 102 2.8(20)) 

-^- 2 F 1 (a,(3; T ,z) = ^ 2 F x {a + 1, + 1; 7 + 1; z). (43) 

az 7 

and the following elementary relation holds (see [12] vol.1 p. 105 2.9(2)) 

2 F l (a,P; r ,z) = (1 - zy- a -P 2 F 1 ( 1 - a, 7 - for,*)- (44) 
The Legendre functions are solutions of Legendre's differential equation 

(1 - z 2 )u"(z) - 2zu'(z) + [a(a + 1) - 6 2 (1 - z 2 ) -1 ]^) = °- ( 45 ) 



Making appropriate substitution it can be reduced to the hypergeometric equation (|4ip 
and consequently its solutions are given in terms of hypergeometric function. More pre- 
cisely, the Legendre function of the first and second kind are defined by (see [12] vol.1 
p.122 3.2(3) and p.143 3.4(6)) 



P aW = T^T-^r^y 72 2^-0,0 + 1; 1-6;^), .»■>-!. 1 .10) 



r(i - b) \\i - x\ 

Q b (x) - e ^ 1/2 na + b+1) b/ 

««W ~ 2a+ l xa+b+1 r(a + 3/2 ) & L > X ^ 



,a + b a + b + 1 , 1 . 

X2-Fi(— + 1, 2 ;a + 3/2;^), x > 1, 



respectively. The functions P^ 1 and are independent solutions of 

6.3 Skew-product representation of Brownian motion 

We now introduce the spherical Brownian motion on the unit sphere S 1 ™ -1 C W n as a 
diffusion on 5™ _1 with the generator being the one half of the Laplace-Beltrami operator 
A^n-i of the manifold 5™ _1 . It is well-known that the following formula holds 

Agn-! = (Sin0) 2 "" A[ (sin0 )«-2 9 + ( Sin ^-2 A 

where the angle is the angle between the pole and the given point on the sphere and 
Agi = Now, if we consider the action of A^n-i on function depending only on <fi, 

this reduces to the generator of the Legendre process LEG{d): 

iA sr , = £ K .m>"- = i ^ + H_l co. * A 

Changing variable cos = t we obtain 

1 - i 2 <9 2 n - 1 9 
2 dt 2 2~ t dt' 

We now invoque the classical skew-product representation of the n-dimensional Brownian 
motion (see, e.g. Ito-McKeane 7.15) stating that it can be represented as the product of 
= {r[ ,t ^ 0} - the Bessel process BES(n), v = n/2 — 1 with the generator 

I d 2 n-1 d 
2or z 2r or 
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and independent spherical Brownian motion = {Q(t);t ^ 0} on S™ 1 with time changed 
according to the formula 

[ ) Jo (R^(s)r 

Moreover, for i / 0, we introduce the process 5 = {S(t);t ^ 0} defined by S(t) = 
b| • ^ ne P rocess & describes the cosine of the angle between the starting point x 
and the spherical Brownian motion 0. Consequently, the cosine between the starting 
point x and W(t) is just S(A(t)). The skew-product representation and the previous 
given considerations imply that S is independent from the Bessel process J?M and the 
generator of S is given by 

„ 1 - t 2 d 2 n - 1 d 

with a domain Dg = {u G C 2 [— 1, 1]; u'(—l) = u'(l) = 0}. Three basic characteristics of 
the diffusion: the speed measure, the scale function and the killing measure are described 
by the following relations (see also [B]) m(dx) = 2(l — x 2 )( n ~ 3 ^ 2 dx, s'(x) = (1 — x 2 )^~ n ^ 2 , 
k{dx) = 0. Moreover, the points —1 and 1 are non-singular reflecting points. We denote 
by pf(x,y) the transition density function with respect to the speed measure, i.e. 

P x (S(t)GA)= f pf(x,y)m(dx), AeB[-l,l]. (50) 
J A 
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